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We investigate effects of electron-electron interaction on the transmission probabil- 
ity of electrons through a tunnel junction in a strong magnetic field. We start with the 
Hartree-Fock approximation, and we show that the coulomb interaction, which gives rise 
to the divergence of Fock correction, should be replaced by the dynamically screened 
interaction. We also show that we should use the bare coulomb interaction for Hartree 
term. We take into account higher order contributions using a simple renormalization 
group approach. The temperature dependence of the transmission probability is quali- 
tatively similar to that of a one-dimensional system. 

KEYWORDS: Friedel oscillations, strong magnetic field, one-dimensional system, screened coulomb inter- 
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1. Introduction 

Many-body effects on the electron transport in one-dimensional systems have been studied by 
several authors. 1-4 ) The transport of interacting spinless Fermions through a single barrier was 
studied by Kane and Fisher 1 ) and the case of spin-1/2 Fermions was investigated by Furusaki and 
Nagaosa. 2 ) They found that the transmission probability vanishes with a power-law as T — » 0. 
They treated the problem within the framework of the Tomonaga-Luttinger liquid theory, 5 ' 6 ) i.e., 
they neglected the backward scatterings in the electron-electron interaction. 

After that, Matveev et al. developed a different treatment in which the effects of the backward 
scattering can be incorporated. 3, 4 ) They started with the Hartree-Fock approximation and treated 
the logarithmic singularity by a simple renormalization group theory. They found that the temper- 
ature dependence of the transmission probability does not obey a simple power-law. According to 
them, Friedel oscillations of the electron density induced by the barrier potential give an essential 
effect on the electrical conduction in one-dimensional systems. 

Those theories give qualitatively the same results, as for the vanishing transmission probabilities 
at zero-temperature. Experimentally, this tendency was observed in quantum wires. 7 ) In this 
experiment, however, electrons are probably scattered by the disorder such as the fluctuation of 
the width of the wire, and the situation is not as simple as the case of a single barrier. 8,9 ) In fact, in 
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spite of the development of the technique, it is still hard to verify the theoretical predictions using 
artificial one-dimensional systems, because the localization effects also contribute to the decrease 
of the conductance at low temperatures. 

In this paper, we investigate the transport of electrons through a tunnel junction with a strong 
magnetic field perpendicular to it (see Fig. 1): The system is one-dimensional like when only the 
lowest Landau levels is occupied, and we can expect a similar effect. In fact, recently, Biagini et 
al. and Tsai et al. treated this problem and have shown that the transmission probability behaves 
like that of one-dimensional systems. 10,11 ) 

In this system, measurement of the transmission probability is much easier than in one- 
dimensional systems. Moreover, one may be able to extract the interaction effects from the data 
by making use of the magnetic field dependence of the transmission probability. 



B 




Fig. 1. Tunnel junctions in a strong magnetic field. Magnetic field B — (0, 0, B) is perpendicular to the insulator 
thin film. 

In order to do so, it is important to estimate the parameters which determine the temperature 
dependence of the transmission probability in approximations as good as possible. The important 
parameter are V(0), and V(2kp), where V(q) is the Fourier transform of the interaction potential 
V(x) and hp is the Fermi wave number. For coulomb interaction, V(0) is divergent, and in Refs. 
10 and 11 the authors used the static screened coulomb interaction for V(0), and V(2kp). This 
replacement is not appropriate, however, as will be shown in the following parts of this paper. We 
will see that we should use dynamical screened interaction for V(0) and bare coulomb interaction 
for V(2k F ). 

In section 2 we present the calculations of Hartree-Fock correction to the transmission am- 
plitude, and in section 3 we will show the way to take into account the higher order terms using 
renormalization group method. 
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2.1 Non-interacting electrons 

First we investigate the transmission of a noninteracting 3D electron through tunnel junctions 
under a strong magnetic field. We ignore the spin degrees of freedom, for we assume that the 
electrons are fully polarized. We choose the z-axis of the coordinate system along the magnetic 
field, and use Landau gauge A(x) = (0, Bx, 0) for the vector potential. Then the single-electron 
Hamiltonian is of the form 

H = ^-(p + eA(x)) 2 + U(z), (2.1) 

where U(z) is the barrier potential of the insulator thin film, and we assume that the potential 
barrier is localized around z = , i.e., U(z) = for \z\ > a. 

We consider the case when the electrons occupies only the lowest Landau level. It is realized 
if the magnetic field is strong enough so that 

B > - e (2^n 2 e ) l l\ (2.2) 

where n e is the electron density. Then, we can label the eigenstates of the Hamiltonian Hq by a 
two-dimensional vector k = (k y , k z ), and the energy eigenvalues and the wave functions <p^(x) 
are of the forms 

1 h 2 k 2 

4 = -**+-^, (2.3a) 
<pUx)=4 v (x,yK z (z), (2.3b) 



where 



(X + kyi B 2 ) 2 

2£ B 2 



e ik y\ (2.4a) 



{ t e tk * z , {z>a), 



i e iK z z + rQe ik z z^ (z > a), 



Here oj c = eB/m, £b = yh/eB, and to an d r o ar e the transmission and reflection amplitudes 
through the barrier, respectively. We assume an infinite system and neglect the effects of the 
boundaries. 
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2.2 First Born approximation 

Next we calculate the correction to the transmission probability due to the electron-electron 
interaction. We calculate it to the lowest order in the interaction within the Hartree-Fock theory. 
Using the electron field operator ip(x), we can write the many-body Hamiltonian H as 



with 



TCo = J dxip\x) 



7~t — Tio + 7ii, 



— (p + eA(x)) 2 + U(z) 



(2.5) 



ifj(x), 



(2.6) 



and 



Hi = ^ J J dxdx'^{x)^{x')V{x - x')^{x')^{x), 



(2.7) 



V{x) being the coulomb interaction potential, and if we do not write explicitly the ranges of the 
space integrations they mean integrations over the whole space region. We calculate the correction 
to the wave functions using Green's function method. 12 ) 

The single-electron wave functions (pk(x,r) can be written in the form 



ip k (x,T) = i f dyG R (x,T;y,T )(p k (y,T ), r > r 



(2.8) 



in terms of the retarded Green's function Gr(x, t; y, To). Assuming that the interaction is switched 
on at time To, we put ifk(x, To) = (p® k (x)e~' l ^k T v/ r \ £° being the single-electron energy measured from 
the Fermi level, and later we will let To — > — oo. The retarded Green's function Gn(x,y,u>) can be 
obtained from the Matsubara Green's function Q(x,y,u> n ) in terms of the analytic continuation 



G R (x,y,uj) =G(x,y,-iuj + r]) 



(2.9) 



where oj n = nkBT(2n+l)/h, n being an integer, and n — > +0. The unperturbed Matsubara Green's 

i. 
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Fig. 2. Feynman diagram for the first Born approximation within Hartree-Fock theory: The thin solid lines indicate 
the Green's function without the interaction correction, and the wavy lines indicate the interaction. 
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function Q°(x,y,u} n ) is denned by 

G{x,y tUn ) = J—^—^. (2-10) 

Hereafter, if we do not write explicitly the ranges of the wave number integrations, they mean 
integrations over the whole wave number region. The Feynman diagrams for Q(x,y,tu n ) to the 
lowest order in the interaction is shown in Fig. 2, and the expression corresponding to it is easily 
found to be 

Q(x,y,uj n ) = Q°(x,y,uj n ) + J dx\dx' x 

x Q c \x,x 1 ,u n )T lHF (x 1 ,x' 1 )g c \x[,y,uj n ), (2.11) 

where Hhf(&i,Xi) is the Hartree-Fock self-energy which is independent of the frequency uo n and 
is given by 

hT, HF (xi,x' 1 ) =S(xi - x[) J dx 2 V(x 1 - x 2 ) 

X k B T^2g (X2,X 2 ,LO n >) 

n> 

- V( Xl - xiJfcBT ^g°(a!i,xi,w n 0. (2.12) 

n' 

Then we obtain the following form for single-electron wave functions (see Appendix A) 

- jj dx 1 dx' 1 G k (x;x 1 )V F (x 1 ,x' 1 )^° k (x[), (2.13) 



where 



f dk' Vu>{x)Wu*{xi) 



and 



V H ( Xl ) = J dx 2 V( Xl - x 2 )p(x 2 ), (2.15) 

V F (x 1 ,x' 1 ) = V(x 1 -x[)J^ ^^(a^Vi), (2.16) 

pix) being the electron density and the subscript oc. means the integration over the occupied 
states, i.e., \k z \ < k F . As for the treatment of the singularity in eq. (2.14), the reader is referred 
to Appendix B. 
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To calculate the correction to the transmission probability we need only the asymptotic form 
of the Gk(x,xi) at z — > +00. From eq. (2.14), we obtain (see Appendix B) 



Gfc(x;xi) = 



2'KihVFlB 



exp 



(x - xi) 2 + {y- yi) 



x exp 



i(x + xi){y - yi) 



2£ B 2 

tQe ik z (z- Zl ) ^ Zl <0, 
e ik z (z- Zl ) +rQe ik z (z+z 1 ) j Zl> Q, 



(2.17) 



with vf = hkF/m. The electron density p(x) in eq. (2.15) is calculated from the wave functions 
(2.3b) with (2.4), and at large distances \z\ it behaves like 

dk 



p(x) ee [ 



oc. (2tt) 

1 



2tt£ b 2 

o " n ( z ) 

2ir£ B 2 



l r ol . kp 
—r—r sm(2fcp \z\ + argr ) H 

2-K\z\ IT 



(2.18) 



In the following we will neglect the space independent part in n(z) which gives a constant Hartree 
potential. Such oscillation of the electron density, i.e., the Friedel oscillation, gives an essential 
effect on the transmission probability, as in one-dimensional systems. 3,4 ) 



2.3 Fock correction 

First we calculate the Fock correction ip k (x) to the wave function, i.e., the second term of 
the right side of eq. (2.13) 

fk F \ x ) = ~ J J dx^x'^kix-x^VFix^x'^ip^x^) 

= ~ J J dx 1 dx 1 G ky (x,y;x 1 ,y 1 )G kz (z;z 1 ) 

dxidyidx'^y^Gky (x, y; x 1 ,y 1 ) 




x 



/dk' 



x 1 1 dzxdz^Gk^z; zi)ul z {z[) 
rkF dk 

-kp 

where G ky (x, y; x\, y±) and G kz {z\z\) are defined in Appendix B. 



/k F jy 



*i)«kVi) 



(2-19) 
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As is shown in the Appendix C, the integration over z\, z[, and k' z gives 

dzidz^Gk^zi) 



//' 



/ k F rlh' 
k f*«2,(*0«2i*(*i) 



t \r \ 2 
2irhvp 



V(xi - x[,y! - yi;0)ln 



1 



\k z - k F \d 



r,ik z z 



(2.20) 



where V(x\ — x[, y\ — y[; q z ) is the Fourier transform of V{x\ — x[) with respect only to z component 
of xi — x[ and is defined by, 



V(xi - x'^yi - y[;k z ) = 

and d is the cut-off length. 

From Eqs. (2.19) and (2.20), we have 



dz' 2 V( Xl - xi, yi - y[, z' 2 )e- ik * z * 



(2.21) 



<Ph '{x)=- 




dxidyidx^dy^G^ (x, y; x 1 ,y 1 ] 



dh' 

— <^(xi,yi)<^ (x 1 ,y 1 ) 
x t J^V s (x 1 -x' 1 ,y 1 -y[;0) 



l0 (J J 



2-Khvp 



x In 



1 



Ak z z 



\k z - k F \d / 

Thus the contribution from the Fock term <^ 7 (a?) is written as (see Appendix C) 
where 



(2.22) 



02(5) = 



1 



?±dgj. e -^ fl2 / 2 V r ( 9± ), 



(2.23) 



(2.24) 



2n 2 hvF Jo 

V(qj_) being the Fourier transform of interaction V{x) with q_|_ = (q x ,Qy,ty, and we obtain the 
following lst-order correction to the transmission amplitude by the Fock term 



t ( 1F ) = -a 2 (B)t (l-\t \ 2 )ln 



\k z — k F \d J 

It is to be noted that if V(x) is coulomb interaction we have 



V(q) 



(2.25) 



(2.26) 



e being the dielectric permittivity of the matter, and that the right hand side of eq. (2.24) is 
divergent. In this context, one of the authors (A.K.) reinvestigated the problem from a point of 
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view different from those of the theories mentioned above: 13 ) Using the parquet diagram method 
developed by Bychkov et a/, 14 ) he has shown that, to be consistent, this singularity should be 
renormalized, too, and that the interaction potential should be replaced by the screened interaction 
within Random Phase Approximation (RPA). 15 ) 

As the screened interaction is frequency dependent, the second term of the right hand side of 
eq. (2.12) should be replaced by 



y~V(a:i - x[,iu n - w nl )Q°(x' 1 ,x 1 ;ij nl ) . 



(2.27) 



According to the arguments in Ref. 13, relevant energy is frvFq z ~ fi^n ~ k B T and we have to 
replace V{q A _) in eq. (2.24) with 

lim V s {q,u n = v F q z ) • (2.28) 
where frequency dependent screened interaction is given by eq. (D-15). Then we easily find that 

-qH B */2 



poo 

a 2 (B)=2(K£ B ) 2 / dq- 



qe 



(2.29) 



■ q 2 + K 2 e -qH B */2 ' 

k being given by eq. (D-18). We easily find that this k is smaller than that in Ref. 11 by a factor 
1/V2. 



2.4 Hartree correction 

Next we calculate the correction to the transmission amplitude by the Hartree term. We can 
write the contribution of the Hartree term ip k (x) as 

( P k 1H \ x )= J dxiG fc (x;xi)V^(a5i)<Pfc(xi) 

dx 1 G ky (x,y;x 1 ,y 1 )G kz (z;z 1 )(()l y (x 1 ,y 1 ) 



J dx 2 V(x 1 - x 2 )p(x 2 ) 



4M) 



= J J dxidyiG ky (x, y; x u yi)4>° ky Vi) 

dx 2 dy 2 f 

dziG kz (z;zi)u kz (zi) 



2tt£ b 2 
j dz 2 V(xi - x 2 )n(z 2 ) 



(2.30) 



Here we should not replace V(x) with the screened interaction according to the following reason. 

The above mentioned replacement can be expressed by Feynman diagrams in Fig. 3 as the 
renormalization of the interaction. On the other hand, this term can be interpreted as the correction 
to the one-electron Green's function as is shown in Fig. 4. In the section 3 we take into account 
the higher order correction, which can be interpreted as the corrections to the one-electron states. 
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Fig. 3. The Hartree correction in terms of the screened coulomb interaction 
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Fig. 4. The Hartree correction in terms of the one-electron Green's function 

Therefore, if we replace the coulomb interaction with the screened interaction, it gives rise to the 
overcounting of the corrections. 

In other words, the effects of the Hartree terms are essentially the screening of the barrier 
potential, and the screening should not be counted twice. Thus, in the Hartree term, we should 
keep the bare coulomb potential in eq. (2.12). 

As is shown in the Appendix C, the integration over z\ and z 2 in eq. (2.30) gives 

<-lziG kz (z;zi)ul z (z{) 
dz2V(xi - x 2 )n(z2) 



27rhvF 
x In 



V(xi - £2,2/1 - V2^k F ) 



\k z — kp\d 



(2.31) 



From Eqs. (2.30) and (2.31), we have 



dxidyiG ky (x,y;xi,yi) 



x 



dx 2 dyr 



4> k Axi,yi] 



x to\r \ ^ _ ^ ^ _ ^ ^ 
ZttTivf 
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1 



hence the contribution from the Hartree term ip^ H ^ (x) is written as (see Appendix C) 

= ai ( B )Ml-M 2 ) in (^^) 

x<>,!/)e il -, (2.33) 

where 

°' (B)S 4$- (2 ' 34) 
Thus we obtain the following lst-order correction to the transmission amplitude by the Hartree 

term. 

t <->= Ql(B )Ml-M 2 )ln(^^). (2.35) 

Eqs. (2.25) and (2.35) give the lst-order correction to the transmission amplitude within 
Hartree-Fock theory 

t« = -a(B)t (l - M 2 ) In ( u \ n ) , (2-36) 



\k z - k F \dJ ' 

where dimensionless parameter a(B) of the electron-electron interaction is 

a(B) = a 2 (B)-a 1 (B), (2.37) 
and the lst-order correction to the transmission probability is given by 

Since this is the result obtained by perturbation theory, it is applicable as long as 

a(B)ln ( l —— ) < 1. (2.39) 



\k z - k F \d; 

However eq. (2.39) is no longer valid at low temperature since only electrons of k z ~ k F contribute 
electric conduction at low temperatures. Thus we have to take into account the higher order 
contributions in the interaction. 

3. Higher order contributions 

In order to include the higher order corrections, we use a simple renormalization group (RG) 
approach called the poor man's scaling developed by Anderson for the Kondo problem. 16 ) 

We assume that only electrons in the strip of the wave number k z of halfwidth Ao = l/d near 
the Fermi wave number k F contribute to the correction (2.38) (see Fig. 5). In this strip we linearize 
the dispersion relation of electron energy: 

ti = hv F {\k z \-k F ). (3.1) 
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(l-s)A« (l-s)ho 



-Ao she k z =0 sA, Ao 

Fig. 5. The reduction of the cutoff Ao. k z is the wave number measured from ftp. 

Note that d is the cutoff of integrations not only due to the validity of eq. (2.18) but also due the 
validity of the linearized dispersion. 

Next we reduce the cutoff Ao to sAq with 1 — s -C 1 (see Fig. 5). If we simultaneously 
renormalize the transmission probability Tq so that the effects of the states excluded by this RG 
transformation are taken into account, the new problem by this RG transformation is equivalent 
to the original problem. The change in Tq found in the first Born approximation is (see Appendix 
E) 

ST = -2a(B)T (l-T )]n(^\, (3.2) 

where E = hvps Aq, Eq = HvfAq are respectively the electron energy measured from Fermi level. 

We apply the RG transformation again, reducing the bandwidth E —* sE step by step. Dur- 
ing each step of rescaling the cutoff, transmission probability T is renormalized according to eq. 
(3.2) with Tq being substituted by the modified T from the previous step. The effect of these 
renormalizations may be found as a solution of the differential equation (RG equation) 

WW) = - MB)r(m - T{m <3 ' 3) 

If we integrate this RG equation from E = Eq to E = £ with the initial condition T\e=e = ^Oi 
the transmission probability becomes 

T (W°(*) 

{) Ko + T (£/E )MB)> W 

where TZq = 1 — Tq. At low temperature £ in eq. (3.4) should be replaced by ksT and the 
temperature dependence of the transmission probability is found to be 

TIT) = (3 5) 

1 ' n + T (k B T/E )^)- (6 - 0) 

Our result is of the same form as that of one-dimensional electron systems except for the parameter 
of the electron-electron interaction. 



3.1 Temperature dependence of the transmission probability 

From eq. (3.5) we find that the transmission probability vanishes as T — > if a(B) > 0, as in the 
case of one-dimensional systems. On the other hand, in Fig. 6 and Fig. 7, we see that a(B) can be 
negative for large magnetic fields. In this case, according to eq. (3.5), the transmission probability 
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should become 1 at OK. Since the calculation of a(B) is based on random phase approximation, it 
is not clear whether the values of a{B) are quantitatively reliable. Nevertheless, the tendency that 
a(B) decreases as the magnetic field increases should be reliable. As can be seem from eq. (2.37) 
the Fock term suppresses the transmission while the Hartree term enhances it. The role of the 
Hartree term can be interpreted as the screening of the barrier potential, and it is not surprising if 
the barrier becomes transparent because of the strong screening due to the one-dimensional nature 
of the system. 



a (S) 
i 




Fig. 6. In the case of n = 10 23 m~ 3 , parameter a(B) of interaction vs. magnetic field B. (i) e = \Qeo,m — 0.5m e , 
(ii) e = 11.9eo, m = 0.26m e , (iii) e = 13.1eo, m = 0.067m e (typical values of GaAs), and (iv) e — 15eo,m = 0.05m e , 
where eo is the permittivity of vacuum and m e is the rest mass of electron. 



a (B) 
A 




Fig. 7. In the case of n = 10 24 m 3 , parameter a(B) of interaction vs. magnetic field B. (i) e = lQeo^m — 0.5m e , 
(ii) e = 11.9eo,7n = 0.26m e , (iii) e = 12>.leo,m = 0.067m e (typical values of GaAs), and (iv) e — 15eo,m = 0.05m e . 



4. Summaries 

We investigated the effects of electron-electron interaction on the transmission of electrons 
through a tunnel barrier in a strong magnetic field. We have found that the temperature dependence 
of the transmission probability is at large the same as that of one-dimensional systems, as in Refs. 10 
and 11. In fact, such behavior can be expected because in both cases the Friedel oscillations of the 
electron density play an essential role. 

In estimating the parameter a(B), we have shown that we should use dynamical screened 
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coulomb interaction and bare coulomb interaction for Fock correction and Hartree correction, 
respectively: It is important to use approximations as good as possible, in order to verify the 
interaction effect from the experimental data. 

It is interesting to investigate the magnetic field dependence of a(B). To do so, one need not 
go to very low temperature. In fact, at moderate temperatures from eq. (3.2) we can expect that 



If a{B) decreases with the increase of the magnetic field, we can expect it to become negative for 
a larger magnetic field, as in Figs. 6 and 7. 

Some line in these figures are for the values of effective mass etc. of typical semiconductor. In 
the case of doped semiconductors, the scatterings of electrons by impurities will suppress the long 
tail of the Friedel oscillations. In order to observe clear interaction effects, the mean free path of 
the electrons has to be much longer than the Fermi wave length. 

Very pure semimetals are good candidates for the observation of the interaction effects. In 
semimetals, however, the band structures are generally very complex, and the theory needs some 
modifications to be applied to them. 
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Appendix A: The derivation of eq. (2.13) by Green's function method 

First we consider only the Hartree term. From Eqs. (2.11) and (2.12), the corresponding part 
of the Matsubara Green's function Q H (x,y,u n ) is given by 



p{x2) being the electron density. Thus the corresponding retarded Green's function can be written 

as 



ST = 2a(B)T (l-T )ln{k B T). 



(4.1) 




(A-l) 



where Vh{x\) is the Hartree-potential 




(A-2) 




<A»( x i) i Pk" (g) 



(A-3) 
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The Fourier transform of eq. (A-3) is given by 

Gr(x,V,T- T ) =i J ^2 / 7^)2 / ^l^fc'(^)^fc'*(^l)Vff(^l)^fc»(^l)^fc»*(2/) 

e -MT-T ) 

X J to(u-e k ,/h+ir,)(u-§ l /h + ir,y ( ' } 

The correction to the wave functions ip k (x, r) by the Hartree term is given by 

(p k 1H \x,r) =i J dyG^(x,y,T - T )<p° k (y , t ) , (A-5) 

where </?fc(y,ro) = (p k {y)e~' l ^ T o/ h an d we will let To — ► — oo. Putting eq. (A-4) into the above 
expression and performing the integral over y, k" , and u, we obtain 

Vfc 1H) («»T) =\ j j dxwl^x^l^ix^VHix^tplix^To) 

In the limit to — »■ — oo, the second term in the curly bracket of the above equation oscillates very 
rapidly when the integral over k! is done, hence its contribution can be neglected. Therefore we 
have 

v rw) = «-«"* / / ^ ^^V i-^c,). (a.7) 

Hereafter we omit the time-dependence of the wave functions, and it can be written in the form 

Vk H) ( x ) = Jd Xl G k (x; x^Vffix^Kx!), (A-8) 

where G k (x;xi) is the single-electron Green's function for noninteracting electrons and is defined 

as 

As for the treatment of the singularity of the integrand, the reader is referred to the Appendix B. 

Next we consider the Fock term. Using Eqs. (2.11) and (2.12), we can write the corresponding 
part of the Matsubara Green's function Q F (x,y,u> n ) as 

g F (x,y,uj n ) = -- J J dx 1 dx' i g (x,x 1 ,io n )V F (x 1 ,x' l )g°(x' 1 ,y,uj n ) (A-10) 

where Vf(xi,x' 1 ) is the Fock-potential 

/dk' 
^¥&(*i)AVi). (A-ll) 
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Thus the corresponding retarded Green's function can be written as 
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The Fourier transform of eq. (A- 12) is given by 



V F (x 1 ,x 1 ) 



, ,V k »{xi)y k „* {y) 



(A-12) 



did 



-jaj(r-ro) 



(A-13) 



Thus, as in the case of Hartree term, we can calculate the correction to f^(x, r) by the Fock term 
from 

¥>k F \x,T)=i J dyG^(x,y,r - T )ip k (y,T ) , (A- 14) 

i.e., 

^ F) (x) = - J J dx 1 dx' 1 G k (x;x 1 )V F (x 1 ,x' 1 )cpl(x' 1 ), (A-15) 

where we omit the time-dependence of the wave functions and G k (x;xi) is defined by eq. (A-9). 
From Eqs. (A-8) and (A-15), we obtain eq. (2.13). 



Appendix B: The calculation of the single-electron Green's function 

To calculate the correction to the transmission probability we need only the asymptotic form 
of the G k (x,x\) at z — > +00. From eqs. (2.14), (2.3b) and (2.4) we have 



G k (x,xi 



dk'o, wo *, ^ [ dki z ul, z {z)ul,*{ Zl ) 
-^-4,( x ,y)4, ( Xl , yi )- « 



= G ky (x,y;xi,yi) ■ G kz (z;zi), (B-l) 
Note that the energy eigenvalues depend only on fc 2 and first we calculate G}~ y (x,y;xi,yi): 

/dk' 
-^^(x,y)^/(xi,yi) 



1 , 

/ y_ P ik y(y-yi) 



x exp 



We let L x — > 00, then we obtain 

GkSx^y-^x^yi) 



L X /2£ B 2 2vr 

(x + fc;^ 2 ) 2 + (xx + fc;V) 2 

2^b 2 

(x - xi) 2 + (y- yxf 



(B-2) 



2vr^' 



x exp 



■ exp 



4^ 



i(x + Xl)(j/ - j/i) 

2/,r' 



(B-3) 
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Next we calculate Gk z (z;z±) for k z > 0: 

In treating the singularities of the integral, we replace k z with k z + irj (77 — > +0) so that we extract 
only the out going waves at z — ► +00. Then, for Z\ > 0, from eqs. (2.4b) and (2.4c) we have 

, dk' z \t \ 2 e ik '^ z - z ^ 



2tt hvF(k z + irj — k' z ) 

r0 rffc / ^ik' z z + ^-ik'^-ik'^ + r(j * e iKz^ 

+ i_ 00 ^T hv F (k z + i v + k' z ) ' ! , " ,) 

where we used the linearized energy dispersion relation eq. (3.1). Since the integrands of the above 
integrals oscillate very rapidly for large z, the main contributions come from only such k' z that the 
denominators vanish. Therefore we may extend the regions of the integrals to (—00,00) and, after 
adding a contour of an infinitely large semi-circle, we obtain (see Fig. 8) 

G '^ (2;Zl) = ^{ eifcz( ^ l)+roeifcz( ^ l) } ' (B ' 6) 





" » \ 


- R Ir X r, 





Fig. 8. The contour of the first integral in eq. (B-5) (R — > 00) 



Similarly, for z\ < 0, we obtain 



1 



ihvp 



t e 



ik z (z—z±) 



where we used the relation 

tor * + t *r = 0, 

which can be derived from a time reversal symmetry argument. 

Thus from eqs. (B-l), (B-3), (B-6) and (B-7) we obtain eq. (2.17). 



(B-7) 
(B-8) 



Appendix C: The calculations of the correction to the transmission amplitude by the 
perturbation theory 

C.l The derivation of eq. (2.20) 

In eq. (2.20), because of V{x\ — x± ), main contribution comes from the regions where z\ and 
z[ are of the same sign. First we consider only the integration over k' z in the Fock potential. When 
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both z\ and z[ are positive, it is easily found that 
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^ d ^l,{ Zl )ul,\z[) 
Jo 27T 



+ ko 



'° rib' 
|2 i CIK Z -ik'^zi-z'j) 

2ir 



kp 







dk' 



s'm[kp(zi — z[)] 



kp 



dk' z 
2tT 



TO 



* e ik' z (z 1 +z' 1 ) +cc 



+ ko 



7r(zi - 4) 

sin[/ci?(zi + z^) + argro] — sin(argro) 



tt(z 1 + z[) 

Similarly, when both z\ and z' x are negative, it is easily found that 
r k * dk' 



(C-l) 



-kp 



v z„,0 



«w(«i)«fci*(2:i) 



2vr 



sin[A;ir(zi — z[)] . sm[k F (zi + z[) — argr ] + sin(argr ) 
+ l r o| 



(C-2) 



7r(zi-^) ' ir(zi + z[) 

Next we transform z\ and into Z2 = z\ + zj and z' 2 = z\ — z[, respectively. Then the left side of 
eq. (2.20) becomes 



to r o\ r o\ c ik z z 
2irihvF 

(•OO 

X 



/OO 
dz' 2 V(xi - x[,y 1 - 
-oo 







^okol 

irihvF 



OO 

e ik z Z2 r0 e 

dz2 sm(kFZ2 + argro) + / dz2 

Z2 J-oo z 2 

roo e ik z Z2 



ik z Z2 



■ s'm(kFZ2 — argro) 



? kzZ V(xi - x[, yi - y[; 0) / dz 2 sm(k F z 2 + argr ), 

JO Z2 



(C-3) 



where V{x\ — x[, yi—y[; k z ) is the Fourier transform for only z component and defined by eq. (2.21) 
Since the expression for the electron density, eq. (2.18), is only valid at large distances \z\, we 
restrict the integration in eq. (C-3) in the regions \z\ > d. At low temperatures, the electronic 
conduction is determined entirely by electrons with energy close to the Fermi energy. Therefore in 
the integral in eq. (C-3), we leave only the terms which are divergent for k z — > kp. Then, we find 
that the integral becomes 



dZ2 



cos[(k z - k F )z 2 ] 

Z2 



and that for \k z — kp\d <C 1 eq. (C-3) reduces to 

*oM 2 T ~, 



2irhvF 



^ 1 -, Ul -, i; ) e-ln(^-^). 



(C-4) 



(C-5) 
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C.2 The derivation of eq. (2.23) 

The integral over k' y in eq. (2.22) is easily found to be 



/ 



elk' 



2tt£ b - 



■ exp 



o */•„/ „./ 

y 

M2 i „/\2l 



{xi _ x ' i) 2 + {yi _ yli y 



exp 



Then we find that eq. (2.22) is given by 



4 1F) (-) 



*okol 

2irflVF 



In 



„ik z z 



/ i 

x < exp 

I 2tt£ b 2 



k z -k F \d; JJJJ-OO 

(xi -xi) 2 + (yi -yi^ 2 



ifol +xj)(j/l - j/j) 



dxidyidx' 1 dy[ 

i{xi + ari)(yi - y{) 



1 



2ttV 
1 



exp 



(xi - sir + -yi) 



7TV4V/2 



exp 



(xj + fc,V) 2 
2^b 2 



exp 
exp 



2^ 2 

i{x 1 +x' 1 ){y l - y[) 
2l B 2 



e ^yiV s ( Xl - x[, yi - y' i; 0) , 



and after performing the integrals, we obtain eq. (2.23). 



(C-6) 



(C-7) 



C.3 The derivation of eq. (2.31) 

From eq. (B-6) and (B-7), we have 



/ dziG kz (z;zi) / dz 2 V{x 1 - x 2 )n(z 2 ) 

J — oo J\z 2 \>d 



4M) 



*0 g ifc z z 



ihvp 



/oo /» 
-oo J|z2|>rf 



^0 c »fc zZ 



dz 2 V(xi - x 2 )n(z 2 ) , (C- 



ihl)F \J-oo JO J J\z 2 \>d 

where the integrations over z 2 are done for \z 2 \ > d because the expression for the electron density 
n(z 2 ) is valid only at large distances \z 2 \. In eq. (C-8), we neglect the first term of the right hand 
side because its integrand oscillates and its contribution to the integral is very small. Then the 
right hand side of eq. (C-8) is easily found to be 

^-e lkzZ [ dz 1 cos(2k z z 1 ) + 2i [ dz 1 sm(2k z z 1 ) [ dz 2 V(x 1 - x 2 )n(z 2 ) . (C-9) 
iriVF U-oo Jo J J\z 2 \>d 

Here we consider the integral with cos(2k z zi): 

/oo r 
dzicos(2k z zi) / dz 2 V (xi — x 2 )n{z 2 ) . (CTO) 
-oo J \z 2 \>d 

Changing the integration over z\ to that over z[ = z\ — z 2 and making use of that V(x\ — x 2 ) is 
an even function of the arguments, we easily find that 



h = V(xi - x 2 ,yi - y 2 ; 2k z ) \ dz 2 cos(2k z z 2 )n(z 2 ) , 

J\z 2 \>d 



(C-ll) 
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where we used eq. (2.21). 

In the integral over z 2) we leave only the term which is divergent at k z — > k F . Then from 
cq. (2.18) we easily find that 

1^*0 1 ~ f°° dz 2 

h = --^-V{x 1 -x 2 ,y 1 -y 2 ;2k z ) / — - sin[2(/c 2 - k F )z 2 - argr ] 
2vr Jd z 2 



^V( X1 - x 2 , yi - y 2 ; 2k z ) sin(arg r ) In ^ 



(C-12) 



In the same way, we can calculate the other term in eq. (C-9) and we find that the right side of 
eq. (C-8) becomes 



MgL e * k ' z V{x 1 - x 2 , yi - y 2 ; 2k F ) In 
where we have replaced 2k z with 2k F in the potential. 



1 



\k z - kF\dJ ' 



(C-13) 



C.4 The derivation of eq. (2.33) 

We easily find that V{x\ — x 2 ,y\ — y 2 ; 2k F ) in eq. (2.31) is given by 

2 

V(xi - x 2 ,yi - y2\2k F ) = {(xi - x 2 ) 2 + (yi - y 2 ) 2 } 1/2 ■ 2k F 



(C-14) 



where Kq(t) is the modified Bessel function. Then from eqs. (B-3) and (C-13), we find that eq. 
(2.31) becomes 



{1H) ( , *oko| 2 e 2 . 
fk (x) = — —In 



2ixhv F 2ixe \\k z — kp\d 



1 




dx 1 dy 1 dx 2 dy 2 < 2 exp 
oo I 2tt£ b 



Mb 



-(x - XlY 



x exp 



1 



L ^B 2 
exp 



{y-yi? 



exp 



2i B 



:(x + x 1 )(y - yi) 



1 



1 ' - - 2\2 



2 ^ 2 V- ^ 



o ik yVl 



rK 



{(*i-* 2 ) 2 + ( yi -y 2 ) 2 } 1/2 .2£: F 



2vr£ B 2 

We first do the integrations over x 2 and y 2 , and using the formula 



poo 
JO 



r = 1, 



we have 



{ih), x *oko| 2 e 2 



2irhv F 2vre Vl&z - k F \d J TKyy ,a/ 4k F 2 £ B 2 



4>lSx,y)- 



J>k z z 



ai(B)t (l- |*o| 2 ) In 



|/c 2 - A; F |d 



<,(*,v)e*-*, 



(C-15) 



(C-16) 



(C-17) 
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where 

1 e 2 

ol\{B) = 



2-Khvp 2vre Ak F 2 £ B 

— 2 . (C-18) 



4k F 

Appendix D: The Screened Coulomb Interaction by RPA 

Below we will consider the case of zero-temperature, and we will not take into account the 
effects of the barrier on the screening, for simplicity. Suppose we put an external charge 

Q(x,t) = Q q e i ^ x -^ (D-l) 

in the system, where q is a three-dimensional vector (below k will be a two-dimensional vector in 
y-z plane, and we define k ± q = (k y ± q y , k z ± q z )). Then the electron density induced by this 
external charge is given by 

Q'(x,t) = x(q,uj)Q(x,t). (D-2) 

Here x(q, w ) is the density response function defined by 17 ) 

1 - f°° 

X (q,u) = -V(q) dt(\p- q ,p q (t)]) (D-3) 
111 Jo 

where (• • • ) indicates the thermal average, V(q) is the Fourier component of the coulomb interaction 
potential 

V(q) = ^, (D-4) 
e being the dielectric permittivity of the matter, and p q is the electron density operator 

Pq = j dxft (x)4>(x)e- iq - X 

= ^afc-qWexp 

k 

We define the Matsubara Green's function T>(q,uj n ), 

V(q,ij n )= dT( Pq (T)p_ q )e^ T , (D-6) 
J o 

where uj n = 7r/c£T(2n + l)/h, n being an integer. Then x(9> w ) can be expressed in terms of 
V(q,uj n ) by analytic continuation 

x (q, u ) = X!glv(q,-iu + Ti), (D-7) 
and the dielectric function is defined as 

eQ(x,t) e 



' i 2 £ 2 

*-§-9x(2*t, - g„) - -f-(fe 2 + q y 2 ) 



(D-5) 



Many-Body Effects on Transmission through a Tunnel Junction 

From eq. (D-5), we find 

o 2 

— ^(Qx +q y ) 



V{q,u n ) = exp 

22 dr ^a fcl _q t (r)a fcl (r)a fc2+(7 t a fc2 ^ e lw " r exp [i£ s 2 g x (fcij, - ^ - q y ) 



fci,fc 2 
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(D-9) 





Fig. 9. Feynman diagram for w„). 

In calculating T>(q, u n ) using Feynman diagram, among the various terms we retain only those 
corresponding to the diagrams shown in Fig. 9. Then we obtain 

V°{q,u: n ) 



l-V(q)V°(q,u n ) 



(D-10) 



where 



V°(q,u; n ) = exp 



i 2 

~-{<ix 2 + q y 2 ) 



(D-ll) 



with f n = irkBT(2n + l)/ft, n being an integer and J^°(k,u n ) is the Matsubara Green's function 
for free electrons 



It is easy to calculate V°(q,uj n ): 



where /p(£p) is the Fermi distribution function and = q 2 + q 2 . 



(D-12) 



(D-13) 



From eqs. (D-7), (D-8) and (D-10), the analytic continuation of the dielectric function is given 



by 



e(q,LO n )=e\l-V{q)V°(q,LO n ) 
The screened interaction is given by 

V s {q,u> n ) 



q 2 e(q,uj n ) 



(D-14) 



(D-15) 
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As is mentioned in the last part of subsection 2.3, the relevant parameter to the Fock correction 
is given by eq. (2.28). For small q z , from eq. (D-13) we have 



V°(a u ) ~ e ~ lW±/2 2{VFqz? 



(D-16) 



hence we find that 

e 2 

lim V a {q,v F q z ) = p^— _ , (D-17) 



where = (q x , q y , 0) and 



e 2 



\j A7r 2 ehv F e B 2 ' ^ D ' 18 ^ 

Appendix E: The derivation of eq. (3.2) 

First we consider only the Hartree term. The lst-order Hartree correction to the transmitted 
wave with cutoff Ao can be written as 

<p£ H \x,A ) = J dxiG k (x;xi)V H (xi)(p%(xi) 

= jj dx 1 dy 1 G ky (x,y;x 1 ,y 1 )(p ky (x 1 ,y 1 ) J J ^^- u k f (x' , A ) , (E-l) 
where x' = {x\ — x 2 , yi — y2, z) and 



itfcf \x' , A ) = J dziG kz {z;zi) j dz 2 V s {x\ - x 2 )n(z 2 ) 



(e-2) 



Here only uj c H (x f , Aq) is dependent on the cutoff. Within the strip shown in Fig. 5, using linearized 
energy dispersion, the single-electron Green's function for noninteracting electrons can be written 
as (see eq. (B-4)) 

^+ A « dk'uUz)ul\ Zl ) 



G ^Z 1 ) = — - - <T) 

Jk F -A Z7r ?fe ? fe . 



1 r K \ y A^ z)u \^ {zi) 



/—sn 
-Ao 



2irhv F J_ Aq : k z -k' z - k F 
2irhv F J_ sAo z k z -k' z -k F 



— f 



2Txhv F J sAo k z -k' z -k F 
ee G[ A \z- Zl ) + Gf\z- Zl ) + Gf\z- Zl ), (E3) 
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where (A), (B) and (G) denote the regions of — Ao < k' z < — sAq, — sAq < k' z < sAq and sAo < 
k' z < Ao, respectively. Similarly, the electron density n(z 2 ) can be written as 

rkF dk z 



Jk 

-L 

= n (A \z 2 )+n (B \z 2 ). 



_2± L.O 

- SAa dk 



dk 



(E-4) 



Thus Uj, (a;', Aq) is of the form 



u 



(IH) 



(a;', Ao) = y dziG^^z;^) ^ dz 2 ^(aJi - * 2 )n (A) (z 2 ) 



+ y tfeiG^fozi) J dz 2 V{x 1 -x 2 )n^ B \z 2 ) 

+ y dziG|j.f(z;zi) J dz 2 V{ Xl -x 2 )n^ A \z 2 ) 

+ y ds^f (z;*i) j dz 2 V{x l -x 2 )rS B \z 2 ) 

+ y dziG]^(z;zi) y dz 2 ^(a;i-^2)n (A) (z2) 

+ y dziG^^zi) y ^ 2 F(o;i-a; 2 )n( B )(z 2 ) 



«2.(*i) 
«2» 
«2» 

«2» 
«2» 
«2.(*i)- 



(E-5) 



We find that the 4th term of the right hand side in eq. (E-5) can be written as u^ H \xq, sAq), 
hence 

4 m) (*', Ao) = 4 m) (sb', ,Ao) + 8ul H) (*', Ao), 



(E-6) 



where 



8ul H \x',A )^J dz 1 G < £\z;z 1 ) 



dz 2 V s (x x - x 2 )n (A \z 2 ) 
+ y dziG^^z;^) y dz 2 y(a;i-a; 2 )n( B )(2 2 ) 
+ y dziG^f (z;zi) y dz 2 y(a;i-a; 2 )n( A )(z 2 ) 
+ y dziG^^;^) y d« 2 ^(a!i-JB2)n^(z2) 
+ y efeiG|^(z;zi) y d^^(asi - a; 2 )n^(z 2 ) 



«2» 



«2» 
«2» 
«2.(*i) 
«2,(*i)> 



(E-7) 



The 1st and 4th terms of the right hand side in eq. (E-7) are the order of (1 — s) 2 , and we 
neglect them. Here we consider only most divergent terms at low temperature which we used in 
the calculations of the perturbation theory. Then we have 



*4fV,Ao) = t ° i lJ^ 2) V(x 1 -x 2 ,y 1 -y 2 -,2k F )ln 



G) 



„ik z z 



(E-8) 



24 Toshihiro Kubo and Arisato Kawabata 

We put it into eq. (E-l) and perform the integrations of x and y directions. Then the trans- 
mitted wave within 1st Born approximation by Hartree term is written as 

<p™{x) =t e ik ^4 y (x,y) + ^ H \x,A ) 
=t e ik * z <j>° ky (x,y) 

+ ^ H \x,sA ) + ai(B)t„(l - |to| 2 )ln (£) ■ e^ z 4 y (x,y) 

= {to + 5t^}e^4 y (x,y) + ^ H \x,sA ), (E-9) 

where 

6t^ = ai (B)t (l - \t \ 2 ) In Q . (E-10) 
Similarly, the transmitted wave within 1st Born approximation by Fock term is written as 

<pW{x) =t e^ z 4 y (x,y)+^ F) (x,A ) 
=t e^ z 4 y (x,y) 

+ <p£ F \x,sA ) - a 2 (B)t (l - |t | 2 )ln Q ■ e ik * z 4>° ky (x,y) 

ee {t + <5t( 1J? )} e^°>,y) + ^ 1F) (x, a Ao). (E-ll) 

with 

5t^ = -a 2 (B)t (l - \t \ 2 ) In Q . (E-12) 
Therefore the transmitted wave within the present approximation is of the form 

<Pk(x) ={t + 5t}e^ z 4 y (x,y) + (^(a, S A ) + ^(a;, aAo)) , (E-13) 

where 

St = -a(B)t (l - \t \ 2 ) In (-) . (E-14) 



Thus we find that ifk(x) does not change if we reduce cutoff from Ao to sAo replacing to with to + 5t 
at the same time. The transmission probability changes: 

% = \t \ 2 ^% + ST =\t + 5t\ 2 , (E-15) 

and we obtain eq. (3.2), where in stead of Aq, we use the energy Eq = HvfAq,E = HvfsAo. 
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